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S: 

' O ■ Abstract 
^ ■ 

Q ■ We present results of a density functional theory study of MnO in the wurtzite 

^ ' structure. Our motivation is provided by recent experiments reporting ferromagnetism 

in Mn-doped wurtzite-structure ZnO. We find that wurtzite MnO a) is not strongly 
energetically disfavored compared with the ground state rocksalt MnO, b) shows strong 
^ magneto-structural coupling and c) has a piezoelectric response that is larger than that 

, of ZnO. These predictions augur well for the creation of ferromagnetic piezoelectric 

' semiconductors based on Mn-doped ZnO. 



1 Introduction 

m . 

■ The wurtzite structure adopted by ZnO does not have a center of inversion symmetry and 

therefore permits the occurrence of piezoelectricity. Indeed previous studies have revealed 
that the ZnO piezoelectric coefficients are substantial jj because of incomplete cancellation 
between opposing ionic and electronic contributions [2 E]. ZnO is already widely used 
. in varistors for surge protection and as filters for incoming television signals. The large 

I piezoelectric response might provide improved electronic device performance, for example, 

O ■ enhanced and tunable optical devices and optical wave resonators. 

Additional interest in ZnO-based systems has recently been stimulated by predictions 
of ferromagnetism in transition metal-doped ZnO using a simple mean field model 0) and 
also by density functional methods Numerous experimental reports of ferromagnetism 

[3 IHl El El El including some above room temperature 23 have also appeared. Although 
' more recent experimental ^21 and theoretical studies call these earlier observations into 

question, there is a clear need for further research in this field. 

In this paper we investigate the fundamental question of what happens to the piezoelec- 
tric response when Mn is substituted for Zn in the wurtzite structure. The answer is not 
obvious; some wurtzite materials, such as BeO, do not show a large piezoelectric response 
even though it is allowed by the non-centrosymmetricity of the crystal structure. The sit- 
uation is even more complicated in the case of transition-metal substituents because of the 
availability of partially-filled 3d states which can become involved in chemical bonding. Such 
states have been shown to strongly suppress the occurrence of ferroelectricity but their 
influence on piezoelectricity has not been previously investigated. We choose Mn as our 
model transition metal dopant because the 3d band of the Mn^+ ion is exactly half-filled, 
with a gap between the up-spin occupied states and empty down-spin states. Therefore we 
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anticipate a more straightforward band structure than in the case of other transition metals 
in which one of the bands is usually partially filled. Also, such an arrangement is more likely 
to result in an insulating state, which is a prerequisite for the occurrence of piezoelectricity. 

Our primary finding is that wurtzite MnO in fact has a larger piezoelectric response than 
ZnO by upto around 3 times. This remarkable result suggests that, in addition to being 
ferromagnetic, transition metal doped ZnO should have a strong piezoelectric response. De- 
vice applications that have been suggested for such magnetic piezoelectrics, include variable 
transducers with magnetically-tunable piezoelectricity, and electric field or stress-controlled 
ferromagnetic resonance devices 16 . Also, the ability to couple with either the magnetic or 
the electric polarization will offer an extra degree of freedom in the design of conventional 
sensors, actuators, transducers and storage devices. The ability to transform changes in the 
magnetic field into electrical voltage and vice- versa could be used as a magnetic field sensor 
for magnetic field measurement, or in electric current measurement. 

The remainder of this paper is organized as follows. In section II we summarize previous 
calculations on MnO in the usual rocksalt structure, and compare our own results with those 
of the literature. In Sec. Ill we describe the electronic structure of MnO in the wurtzite 
structure. In Sec. IV we present the results of our studies of magneto-structural coupling 
and piezoelectricity in wurtzite MnO. In the last section we present our conclusions and 
suggestions for future work. 

2 MnO in the Rocksalt structure 

The experimental ground state crystal structure of MnO is rocksalt \17\ . In this section we 
present results of our calculations for rocksalt-structurc MnO in three magnetic orderings, 
namely the ferromagnetic (FM) and two types of anti-ferromagnetic (AFM) orderings. In 
type-I AFM ordering (AFI), the Mn ions are ferromagnetically aligned in the [001] planes, 
with adjacent planes aligned antiferromagnetically to each other, while in type-II order- 
ing (AFII) they are ferromagnetically aligned in [111] planes, again with antiferromagnetic 
alignment between the planes. The type-II structure is the experimentally observed ground 
state. The purpose of this section is two-fold. First to compare with existing data in order 
to verify our methodology and pseudopotentials, and second as a point of comparison for 
our later studies on MnO in the wurtzite structure. 

Our total energy calculations were carried out using the local spin density approximation 
(LSD A) to density functional theory (DFT) with pseudopotentials and a planewave 

basis as implemented in the ABINIT 20; computer code. We used the optimized pseudopo- 
tentials developed by Rappe et ai, |2] for both the Mn and the O atoms. Our semi-core 
Mn pseudopotential was constructed from a 3s^, 3p^, Mn^+ valence configuration, with 
cutoff radii for the s, p and d states of 1.2, 1.1 and 1.4 a.u.s respectively. Optimization of the 
energy cutoff was carried out using a wave vector, Qc, of 9 a.u.s. with four Bessel functions. 
The oxygen pseudopotential used a 2s^, 2p^ neutral valence configuration with cutoff radii of 
1.2 and 1.4 a.u. for s and p states respectively. Optimization of kinetic energy cutoff for the 
oxygen states used QcS of 7.0 and 6.5 a.u.s with four and three Bessel functions respectively. 
Good convergence was obtained with a plane-wave energy cutoff of 82 Ry for these pseu- 
dopotentials. Rhombohedral magnetic unit cells, containing two formula units of MnO, were 
used for all three magnetic configurations. Highly converged special k-point sets consisting 
of 216 points in the rhombohedral Brillouin zone were used in each case. The exchange- 
correlation functional was parameterized using the Perdew-Zunger parameterization j22] of 
the Ceperley- Alder exchange-correlation functional. 

First we optimized the structures for each magnetic ordering and calculated the total 
energies, which are listed in Tabled As expected, we found the AFII configuration to be 
lowest in energy, with the AFI 0.55 eV/formula unit higher, and the FM 0.60 eV/formula 
unit above AFI. This ordering is consistent with experiment ^T], first-principles Hartree 
Fock calculations [22 and previous LSDA density functional calculations |25[ I26L I77| . Our 
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Magnetic ordering 


Energy (eV/formula unit) 


AFII 


0.00 


AFl 


0.55 


FM 


0.61 



Table 1: Relative energies of different magnetic configurations in rocksalt-structure MnO 



optimized rhombohedral lattice constant for AFII MnO was 8.16 a.u., consistent with pre- 
vious LSDA calculations |2HI and showing the typical LSDA underestimation relative to the 
experimental value of 8.38 a.u. |17| . 

In Figure n we show our calculated band structure of AFII MnO which is calculated at 
the experimental value of 8.38 a.u. (The zero of energy is set to coincide with the highest 
occupied orbital). The occupied d states lie at the top of the valence band, and their splitting 
by the rhombohedral crystal field into two doublets and a singlet can be seen clearly. Below 
and slightly separated from them are the predominantly oxygen p states which have a band 
width of around 4.5 eV. We will show later that the bandwidths of the Mn and O bands are 
affected by both the structure and by the magnetic ordering of the system. We obtain an 
energy gap of 0.75 eV, in perfect agreement with earlier LDA calculations |^. Analysis of 
the orbital character of the LSDA band gap [211 indicates that it lies between predominantly 
up-spin d states at the top of the valence band, and down-spin d states on the same Mn 
atom at the bottom of the conduction band. 




Figure 1: Band structure of AFII MnO plotted along the high symmetry lines of the rhom- 
bohedral Brillouin zone. 



3 MnO in the wurtzite structure 

Having established that our pscudopotentials accurately reproduce previous results for rocksalt- 
structure MnO, we now turn to the main part of our paper, the study of MnO in the wurtzite 
structure. As mentioned in the introduction, our motivation is provided by recent experi- 
ments in which transition metals are doped into ZnO, where their local environment is the 
four-fold coordinated structure considered here. 

We study the six magnetic orderings of wurtzite MnO shown in Figure 13 (Note that 
these were previously studied for MnS by Hines et al. ^Q.) Identical unit cells, containing 
8 formula units, were adopted for all structures, and the pseudopotentials and cut-offs were 
the same as those described in Section |5] A 4 x 4 x 4 k-point grid was used in all the 
calculations. A two-dimensional projection of the unit cell, looking down the c axis is shown 
in Figure 13 The filled circles represent one set of ions in the the z=0 and the z=l planes , 
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while the empty circles show the positions of the other ions in the ■l— i plane. 
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AF2b 
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AF2a 



Figure 2: Different magnetic orderings for wurtzite structure MnO studied in this work. 




z=0,l 
z=l/2 



Figure 3: Two-dimensional projection of the unit cell of wurtzite MnO used in this work 



4 Energetics 

Initially we set the lattice parameters to those of ZnO, and found that the AF3 structure 
was lowest in energy. However prohibitively large stresses were obtained for all the magnetic 
orderings at these lattice parameters. Therefore we chose new lattice parameters, a and 
c of 3.2847 A and 5.1789 A( around 1 % larger than the corresponding values in ZnO) 
and a u value of 0.3687. These new values gave acceptably small stresses. In order to 
decouple the effects of changes in structure from effects of changes in magnetic ordering, 
we then calculated the properties of all our magnetic structures at these lattice parameters, 
and obtained the relative total energies for the wurtzite structures listed in Table |3 The 
calculated energy ordering is 

AF3 < AFl < AF2b < AF4 < AF2a < FM. 

Note that the energies of the AF3 and AFl phases are very similar, as are those of the 
AF2b, AF4 and AF2a. We relate this below to the ferro- and antiferromagnetic interactions 
between Mn nearest neighbors. As in the rocksalt case, the ferromagnetic phase is highest 
in energy. It is notable that the relative energies per formula unit of the AF3 wurtzite 
structure and the AFII rocksalt structure are the same within the errors of the calculations. 
This suggests that it is not unfavorable to incorporate Mn into a tetrahedral environment. 

In Figure 0] we show the band structures for AF3 and FM wurtzite structure MnO (only 
the up-spin states are shown in the FM case). Note that in both cases we obtain metallic 
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Magnetic ordering 


Energy (eV/formula unit) 


AF3 


0.000 


AFl 


0.002 


AF2b 


0.098 


AF4 


0.113 


AF2a 


0.118 


FM 


0.426 



Table 2: Relative energies of the different magnetic orderings in wurtzite MnO. The energies 
are given relative to the reference lowest energy AF3 structure. 



AF3 FM 




rA HK rrA hk r 

Figure 4: Band structures of wurtzite MnO for AFM (left) and FM (right) orderings, only 
the up spin states are shown in the FM case. 

band structures, with the FM phase having numerous Mn d derived bands crossing the Fermi 
energy. The AF3 phase (and the other AFM phases, which are not shown) is close to being 
gapped, with only the Mn 4s band crossing the Fermi energy at the F point. We point out 
that the metallicity is an artifact of the underestimation of the gap by the LSDA, and in 
practice we anticipate that all of the AFM structures should be insulating. Again the top 
of the valence band is composed of predominantly Mn 3d states, which this time show the 
e - t2 crystal field splitting of the tetrahedral coordination. The predominantly oxygen 2p 
bands are again below the Mn d bands and separated from them by a gap. The bottom 
of the conduction band is largely composed of Mn 3d states, except for the broad 4s band 
that crosses the Fermi energy. Compared with the rocksalt band structure, the bands are 
narrower, as anticipated for Mn in four-fold rather than six-fold coordination. In addition, it 
is notable that the wurtzite AFM structure has broader bands than the wurtzite FM. This 
reflects the increase in hybridization possible, via the usual superexchange interaction, in 
the AFM bonded configuration. Mn-0 hybridization is unfavorable in the FM configuration 
between d^ Mn^+ ions as a result of Hund's rule. Cartoons of the Mn-0 bonding in the 
AFM and FM cases are shown in Figure El 
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Figure 5: Schematic of (a) favorable Mn-0 hybridization in the AFM configuration and 
(b) unfavorable hybridization in the FM configuration. In (a) the Mn ions are antiferro- 
magnetically coupled. Therefore the down-spin oxygen 2p electron undergoes covalent bond 
formation with an up-spin Mn 3d electron (left) leaving an up-spin oxygen 2p electron to 
form a covalent bond with a down-spin Mn 3d electron (right). In (b) the Mn ions are 
ferromagnetically coupled. Therefore if the down-spin oxygen 2p electron bonds with an 
up-spin Mn 3d electron on the left (as in the antiferromagnetically coupled case above) the 
remaining up-spin oxygen 2p electron is unable to bond with the up-spin Mn 3d electrons 
on the right. Therefore the amount of hybridization, and hence the band width, is reduced 
in the ferromagnetic case. 



4.1 Exchange Couplings 

To understand the relative energies of the different spin orderings, we next analyze the 
energy differences between the configurations using a two-parameter Heisenberg model; 

H ^ Eo + Y,J^jS,.S, (1) 

where Si.Sj =+1 if the spins i and j are parallel, and -1 if they are antiparallel. We include 
two different values for J; an in-plane Mn-Mn interaction, Jp, and an out-of-plane (along c) 
interaction, Jc, as shown in Figure|B| Note that, for the ideal wurtzite structure, the in-plane 
and out-of-plane Mn-Mn distances are identical, and we would expect to obtain Jp — Jc 
in the absence of any kind of spin-orbit coupling. However, we constrain the parameter u 
to be 0.368 which is smaller than the ideal value of 0.375 and hence we find a percentage 
difference in Jp and Jc of the order of the deviation from the ideal value. 




Jp 



Figure 6: Schematic of the exchange parameters Jp and Jc- 



TableOlgives the number of parallel and antiparallel in-plane and out-of-plane nearest neigh- 
bors for a Mn ion in each magnetic structure. Note that the AF3 and AFl structures have 



6 







AF3 


AFl 


AF2b 


FM 


AF4 


AF2a 


in-planc 


n 


4 


4 








4 


4 




TT 


2 


2 


6 


6 


2 


2 


out-of-plane 


n 


4 


4 


6 





2 


2 




TT 


2 


2 





6 


4 


4 



Table 3: Number of in-plane and out-of-plane FM and AFM interactions for each magnetic 
configuration. 

identical types of nearest-neighbor (nn) Mn-Mn interactions (4 in-plane AFM and 4 out-of- 
plane AFM), which explains their similar energies. Similarly, AF4 and AF2a have identical 
numbers of nn interactions to each other, but have fewer AFM interactions than AF3 and 
AFl (4 in-plane AFM and 2 out-of-plane AFM). AF2b, which is close in energy to AF4 
and AF2a has the same total number of nn AFM interactions (six) , but they all occur out- 
of-plane. Finally the FM phase has no AFM interactions, leading to its correspondingly 
higher energy. By fitting the energies of the six configurations onto the Heisenberg model, 
we obtained the parameters Jp — —0.023 eV and Jc = —0.028 eV, with a root mean square 
deviation of 0.003. The small deviation indicates that a nearest-neighbor Heisenberg model 
describes the energetics of the system well, and suggests that more distant exchange inter- 
actions are weak compared to the nearest-neighbor interactions. Since the Mn-Mn in-plane 
bond length is almost equal to the out-of-plane Mn-Mn bond length, we obtain Jp w Jc as 
expected. 

To conclude this section, we note that AFM coupling is more favourable than FM for both 
in-plane and out-of-plane Mn-O-Mn bonds and. This suggests that, if doping of Mn ions into 
ZnO results in clustering of the magnetic ions, we will obtain strong antiferromagnetism. 

4.2 Forces and stresses 

In many magnetic compounds, a change of magnetic ordering causes a stress which induces 
a structural distortion For example, in rocksalt MnO, AFM spin ordering causes a 

reduction in magnetic symmetry from cubic to rhombohedral, which is accompanied by a 
structural transition j32| . In some materials, in which there is strong magneto-structural 
coupling, the distortion can be quite large. For example in CrN the magnetic ordering 
induces a transition from cubic to orthorhombic symmetry [^OHIISI- And in LaMnOa the 
observed A- type antiferromagnetic ordering is unstable in the presence of distortions |35|. 
The concept of magnetic stress has been introduced to describe structural phase transitions 
that are induced by magnetic ordering |31| . In this section we investigate the magneto- 
structural coupling in wurtzite MnO, as a precursor to studying the piezoelectric properties 
in the next Section. 

Keeping the lattice parameters of all the structures fixed at the AF3 values, the spin 
ordering was changed to each of the other five configurations in turn, and the electronic 
structure was obtained and the stresses and the forces for each magnetic phase were cal- 
culated. In the table below, the relative stresses on the unit cell along the z direction and 
largest forces on the Mn atoms are tabulated for each magnetic ordering. 





AF3 


AFl 


AF2b 


AF4 


AF2a 


Stress <T33 
(Ha/bohr^) 


-7.35 xlO-'^ 


-7.067 xlO-^ 


8.35 xlO-^ 


-2.39 xlO-^ 


-2.46 xlO-4 


Force (max) 
(Ha/bohr) 


5.87 X 10-^ 


5.76 xlO-^ 


4.42 xlO"^ 


8.43 xlO"^ 


8.24 xlO-^ 



We observe similar groupings as those obtained for the relative energies. The AFl and AF3 
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phases are similar in the spin-environment and therefore have similar stresses and forces 
for the same structure. Likewise, the AF2a and AF4 spin orderings show similar stresses 
and forces to each other, but both are considerably larger than those in the AF3 and API 
structures. The significant change in the stress when the magnetic ordering is changed 
indicates a strong coupling between the magnetic moments of the atoms and the crystal 
structure. 

5 Piezoelectricity 

The piezoelectric constant relates the stress response of a material to a uniform electric 
field or conversely the macroscopic polarization induced by a a macroscopic strain. The lat- 
ter definition is used here to calculate the electro-mechanical response from first principles. 
We calculated the changes in macroscopic polarization using the Berry phase approach 
within the ABINIT density functional program 20 . First, the equilibrium structure pa- 
rameters (ao, Co and uq) were obtained for each magnetic configuration, in turn and the 
spontaneous polarizations, Pq, (along the z-axis) were calculated. Then the reference struc- 
ture was strained along the z axis by 1 % (e^z — 0.01, e^z is the percentage change in strain 
along the z-axis ) and its structure relaxed with respect to the internal strain (structural 
parameter u). The polarization Pi for the relaxed structure was then calculated, and the 
piezoelectric co-efficient 633 was obtained using its definition: 

633 = (2) 

^zz 

Subsequent application of a -1% strain along the z-axis showed that the polarization deviated 
slightly from linearity in MnO (in contrast to ZnO, where it is strictly linear over that range 
with the same set of convergence parameters.) All the Berry phase polarization calculations 
used a 2x2x8 k-point mesh and a wavefunction tolerance set to 10~^^. Calculations for one 
of the configurations were repeated with a 3x3x8 mesh to test convergence. A 2x2 mesh 
in the ab-plane is found to be sufficient as it corresponds to a 4x4 mesh for the primitive 
unit cell of the wurtzite structure. In order to avoid the F point, where the gap vanishes 
indicating metallicity, the k-point grid was shifted by ^ along the three co-ordinate axes. 

Our calculated values of 633 for all the five magnetic orderings are listed in Table 3. 
While we find groupings of configurations (AF3 and AFl), (AF4 and AF2a) such that the 
energies and stresses are similar for members of a given group, their piezoelectric constants 
are found to be somewhat dissimilar, as they are second order derivatives of energy. 

The piezoelectric constant of wurtzite MnO in the lowest energy AF3 configuration is 
about 3 times larger than that of ZnO (1.0 — 1.2C/m^|H]). In ZnO, the phonon (internal 
strain) contribution to the piezoelectric constant was found to dominate and this was at- 
tributed partly to hybridization between the filled 3d band of Zn and the 2p band of oxygen. 
In AF3 MnO, we find that the electronic contribution to 633 is also significant and in fact 
has the same sign as that of the phonon contribution (in contrast with ZnO, where the 
two contributions partially compensate). The large electronic contribution arises because 
the d band of Mn is only partially occupied and can therefore contribute directly to the 
piezoelectric response through change in its electron content. 

Note that the enhanced piezoelectric response which we find in wurtzite MnO is con- 
sistent with earlier observations of large dielectric and piezoelectric constants in systems 
that are close to electronic and structural instabilities ^37, ■ Wurtzite structure MnO is both 
unstable structurally (with the rocksalt structure being the grounds state ) and has other 
magnetic configurations close in energy. Finally we point out that Mn doped ZnO should 
have a strong piezoelectric response provided that it remains insulating. 
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Ordering 


AF3 


AFl 


AF2b 


AF4 


AF2a 


633 (C/m^) 


4.81 


11.24 


0.911 


2.05 


1.30 


Po(C/m^) 


0.045 


0.053 


0.017 


0.014 


-0.072 



Table 4: Piezoelectric constant, 633 and the spontaneous polarization, Pq, (C/m^) for the 
different magnetic orderings of wurtzite MnO. 

6 Summary and Conclusions 

In this research we investigated the electronic properties of wurtzite MnO using density 
functional calculations. Using Berry phase polarization calculations we found that mag- 
netic MnO is a stronger piezoelectric than non-magnetic ZnO. This to our knowledge is 
the first calculation of piezoelectric response in a magnetic material and shows clearly that 
presence of magnetism does not preclude a piezoelectric response and can even enhance it. 
The reason for the enhancement is that, in addition to the phonon contribution found in 
ZnO, the electrons in the partially filled d band of Mn also contribute to the piezoelec- 
tricity. Our results suggest that piezoelectricity should persist in Mn doped ZnO provided 
that it remains insulating. On analyzing the energetics of the different magnetic orderings 
for wurtzite MnO, however we find that anti-ferromagnetic ordering is more stable than 
ferromagnetic ordering with the most stable structure having the largest number of nearest 
neighbor antiferromagnetic interactions. From this we can conclude that if doping of Mn in 
ZnO results in clustering, ferromagnetism is unlikely. 
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